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^ ■ We use data on the deep inelastic structure function F2 in order to 

{S| ! constrain the cross-section for scattering a colour dipole off a proton. 

The data seem to prefer parameterisations which include saturation 
effects. That is they indicate that the strong rise with energy of the 
f*"^ , dipole cross-section, which holds for small dipoles, pertains only for 
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Abstract 



r < r s (x) where r s {x) decreases monotonically as x decreases. Subse- 
quent predicitions for the diffractive structure function also hint 
at saturation, although the data are not really sufficiently accurate. 



^ 1 Introduction 

The possible observation of gluon saturation effects in the HERA data at 
very low x-values has been discussed extensively in the context of the colour 
dipole model. In particular, the "saturation model" of Golec-Biernat and 
Wusthoff [H |2] was shown some years ago to give an elegant and accurate 
account of deep inelastic scattering (DIS) at small x, and a rather good de- 
scription of the diffractive deep inelastic scattering (DDIS) data jSJIl]. More 
recently, a new saturation model, the Colour Glass Condensate (CGC) model 
of Iancu, Itakura and Munier jH] , has been formulated. This model, can be 
thought of as a more sophisticated version of the Golec-Biernat-Wusthoff 
model. Gluon saturation effects are now incororated via an approximate 
solution of the Balitsky-Kovchegov equation [S], which applies in the pertur- 
bative region when the gluon densities become large and non-linear effects 
become important. The resulting dipole cross-section, obtained by fitting 
the free parameters of the model to the DIS data at low x, is very similar to 
that obtained in the original saturation model* and Forshaw, Sandapen and 

*For an explicit comparison of the two dipole cross-sections, see jZj. 
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Shaw [S] have shown that it leads to similar, successful, predictions for the 
DDIS data. 

It is clear from this and other work [§] that the predictions of satura- 
tion models are compatible with a wide range of data. However, about the 
same time as the saturation model was formulated, Forshaw, Kerley and 
Shaw (FKS) proposed a two-component Regge dipole model [TU] which ex- 
plicitly excludes gluon saturation effects. As in the saturation model, the 
parameters of the model were tuned to fit the DIS data at low-x, and were 
subsequently shown [TT] to give a good account of the DDIS data without 
further adjustment. Subsequently, the predictions of all three models for the 
DDIS data have been compared to each other jH], showing that they were all 
in good agreement with the data. However, at energies just above the mea- 
sured region, the predictions of the saturation models and the unsaturated 
two component model rapidly diverge. A similar pattern emerges in deeply 
virtual Compton scattering jT^j and vector meson production [7] , although in 
the latter case, there are large ambiguities associated with the wavefunctions 
of the vector mesons. 

From this it is clear that the data analysed to date are insufficient to 
establish the existence of gluon saturation. However, since the saturation 
model and the two-component Regge model were originally formulated, much 
more precise data have become available on both the DIS structure function 
F2 J3| and the DDIS structure function F® 3 [T3]. In this paper we will 
analyse these new data to see if they can throw more light on the question 
of gluon saturation. 

2 Colour Dipoles and Saturation 

In the colour dipole model P^ ll6j. the forward amplitude for virtual Comp- 
ton scattering is assumed to be dominated by the mechanism illustrated in 
Figure ^ in which the photon fluctuates into a qq pair of fixed transverse 
separation r and the quark carries a fraction z of of the incoming photon 
light-cone energy. Using the optical theorem, this leads to 



for the total virtual photon-proton cross- sect ion, where ipL,T(z,r) are the 
appropriate spin-averaged light-cone wavefunctions of the photon and <j(s*, r) 




(1) 
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Figure 1: The colour dipole model for <y*p — > 7*p. 



is the dipole cross-section. The dipole cross-section is usually assumed to be 
independent of z (see below), and is parameterized in terms of an energy 
variable s* which depends on the model. 

In the dipole model, not only the DIS total cross-section, but also the for- 
ward amplitudes for DDIS and DVCS, are determined by the same ingredi- 
ents: the light-cone wavefunctions of the photon and the dipole cross-section. 
We shall briefly discuss each in turn. 

For small r, the light-cone photon wavefunctions are given by the tree 
level QED expressions [11)1 IT?] : 

\Mz,r)\ 2 = - 2 a em Y,e 2 f Q 2 z 2 {l-zfKl{er) (2) 
71 f 

\Mz,r)\ 2 = ^a em J2e}{[z 2 + (l-z) 2 ]e 2 K 2 (er)+m}K 2 (er)} (3) 

where e 2 = z(l — z)Q 2 + m 2 , K and K\ are modified Bessel functions and 
the sum is over quark flavours / with quark masses rrif. These expressions 
decay exponentially at large r, with typical r- values of order Q _1 at large Q 2 
and of order mj 1 at Q 2 = 0. However for large dipoles^ r > lfm, which are 
important at low Q 2 < 4m 2 , a perturbative treatment is not appropriate, 
and hadronic states of this transverse size are more sensibly modelled as 
pairs of constituent quarks, or vector mesons. For this reason, FKS ^Hl 
used constituent quark masses, together with an enhancement factor in the 

^The contributions from very large dipoles r» 1 fm are heavily suppressed at all Q 2 
because of the exponential decay of the wavefunction at large r. 
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region v ~ 1 fm, as suggested by Generalised Vector Dominance (GVD) 
ideas [TBI EH 120] • in contrast, the saturation models retain the perturbative 
wavefunction, but use a much smaller quark mass, which again enhances the 
wavefunction at large r. We have investigated both these approaches and 
find that the difference between them is only important when confronting the 
photoproduction data (mainly that from fixed-target experiments [2JJ ). Since 
we are primarily interested in determining the need for saturation effects 
we will not need to consider these data here, and will therefore use the 
pertubative QED wavefunction (J2J) throughout, treating the quark mass as 
a parameter which can be used to adjust the wavefunction at large r- values. 
At high Q 2 where small dipoles dominate, the wavefunctions are insensitive 
to the quark mass. 

Turning to the dipole cross-section, it is useful to distinguish three re- 
gions: small r, where pertubative ideas are relevant; large r, where typical 
"hadronic" behaviour is expected; and intermediate dipole sizes. The dipole 
formula (JJJ can be derived in the leading log(l/x) approximation of pertur- 
bative QCD [22] • In this approximation, the wavefunction is given by the 
perturbative expression ([2]) and the dipole cross-section is given by 

Att 2 f dk 2 ix 2 r 2 
a(s\r) = —a s j — f(x,k 2 )(l - J (kr)) ~ —a s xg(x, A 2 /r 2 ), (4) 

where f(x, k 2 ) is the unintegrated gluon density and the approximate equal- 
ity holds for sufficiently small r (with A ~ 3). As x gets smaller, the gluon 
distribution grows rapidly, and the dipole cross-section in this region is some- 
times approximated by the "hard pomeron" behaviour 

cr(s*,r) « ar 2 x~ XH (5) 

with Xh ~ 0.4. At large r > 1 fm, which is important as Q 2 — > 0, the models 
we discuss all build in "soft pomeron" behaviour 

a(s*,r)^a s Xs , (6) 

with Xs small or zero. The dipole cross-section is then interpolated between 
the two in the intermediate r region, in a way that depends on the specific 
model. For example, in the original Regge dipole model model of Forshaw, 
Kerley and Shaw [TUl lllj. the dipole cross-section was extracted from DIS 
and real photoabsorption data assuming a form with two terms, each with a 
Regge inspired s dependence: 

a(s, r) = a soft (r)s Xs + a hard (r)s XH (7) 
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where the values Ag ~ 0.06, Xh ~ 0.4 resulting from the fit are characteristic 
of soft and hard diffraction respectively. The functions a so f t (r), a hard {r) were 
chosen so that for small dipoles the hard term dominates yielding a behaviour 
a — > r 2 ( r 2 s) A » as r — > 0. This in accordance with the colour transparency 
ideas embodied in (j3J), since the dimensionless variable r 2 s is closely related 
to x at large Q 2 , where the typical dipole sizes r oc Q" 1 . For large dipoles 
r > 1 fm the soft term dominates with a hadronlike behaviour a « a (r 2 s) Xs . 

In any model which incorporates the "hard pomeron" behaviour (J3J), the 
dipole cross-section increases rapidly with decreasing x for the small dipole 
sizes which dominate at large Q 2 . In nature, however, one might expect 
this sharp rise to be softened by unitarity effects, an effect which we shall 
refer to as "gluon saturation" since, if it occurs in the region where (J3J) is 
valid, it can be interpreted as a dampening of the rapid increase of the gluon 
density with decreasing x implied by (J3J). Such effects are not included in 
the FKS model fOJE]; where for a small dipoles of a given r, the behaviour 

holds indefinitely as x — > 0.* Non-linear saturation dynamics is, however, 
explicitly incorporated into the CGC model, in which the dipole cross-section 
is assumed to be of the form [5] 

„r ln(2/rQ s ) 1 

a(x, r) = 2ttR 2 Mq ( — J for rQ s < 2 



2 

= 2tt R 2 {1 -exp[-a In 2 {brQ s )]} for rQ s > 2 , (8) 

where the saturation scale Q s = (x /x) x ^ 2 GeV and x = Q 2 /(Q 2 + s) is 
the Bjorken scaling variable. The parameters are fixed by a combination of 
theoretical constraints and a fit to DIS data§. This dipole cross-section is 
characterised by a rapid increase with decreasing x, not dissimilar to (0), at 
small r, changing to a softer energy dependence as r increases beyond 

r s (x) = 2/Q s = 2(x/x ) x/2 . 

Saturation arises because of the decrease of the "saturation radius", r s , with 
decreasing x. Let us consider a dipole of fixed transverse size r. If r < r s (x), 
the dipole cross-section increases rapidly as x decreases. However, this rapid 



■fThis is not to deny the inevitable role of nonlinear dynamics, merely to imply that it 
may not be needed in the range of existing data. 

§The authors of [5] carried out a series of fits for different values of the parameter iVo, 
and found that the results were only weakly dependent on this parameter. In this paper 
we choose the fit corresponding to No = 0.7. 
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rise eventually switches to a softer x dependence when x becomes so small 
that r s (x) itself decreases below the the fixed dipole size r. 

Finally we remark that the word "saturation" is also often used to describe 
a somewhat different effect. In all models, the approximate r 2 dependence 
which pertains in the perturbative region at small r is flattened when r 
increases to values in the non-perturbative region. This "non-perturbative 
saturation" is important in describing the transition from high to low Q 2 
values. However it is not what we are concerned with here. 



3 Analysis of the HERA data 

The parameters of the original Regge dipole model of FKS were determined 
using the DIS data available in 1999. However, since then more precise 
measurements of the deep inelastic scattering data in the diffractive region 
have been made at HERA ^3], and there exists new data on DDIS ^3]. In 
this section, we examine these newer data to see whether they can throw any 
new light on the issue of gluon saturation. 

Our strategy is to confront the HERA data with a new Regge dipole 
model (which can be thought of as an update of the original FKS model 
to accommodate the latest data) and also a variant of that model which 
includes saturation. We shall also always compare to the predictions of the 
CGC saturation model. 

Let us introduce our new and very simple Regge inspired dipole model. 
We shall assume that 

a(x m ,r) = Anr 2 x^ H for r < r and 

= A s x^ s for r > n, (9) 

where 

^^r^)' (io) 

For light quark dipoles, the quark mass m is a parameter in the fit, whilst 
for charm quark dipoles the mass is fixed at 1.4 GeV. 

In the intermediate region tq < r < t\ we interpolate linearly between 
the two forms of ( J9*| ). Whether this is a Regge inspired model or a saturation 
model depends entirely upon the way in which the boundary parameter r is 
determined. 
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3.1 Regge Fit to DIS data 



If the boundary parameter tq is kept constant then the parameterisation 
reduces to a sum of two powers, as might be predicted in a two pomeron 
approach. It is plainly unsaturated, with the dipole cross-section obtained at 
small r-values growing rapidly with increasing s at fixed Q 2 (or equivalently 
with decreasing x) without damping of any kind. 

In all cases, we shall fit the electroproduction data in the kinematic range 

0.045 GeV 2 < Q 2 < 45 GeV 2 x < 0.01 , (11) 

which is identical to that which was used to determine the parameters of the 
CGC saturation model. 

The best fit obtained with what we shall henceforth call our FS2004 Regge 
fit is shown as the dashed line in Figure El (left) and the parameter values 
are listed in Table [T] While the values of the Regge exponents Xs and A#, 
and of the boundary parameters ro and r± are eminently sensible, the quality 
of the fit is not good, corresponding to a x 2 /data point of 428/156. This is 
not just a failing of this particular parameterisation. We have attempted to 
fit the data with other Regge inspired models, including the original FKS 
parameterization, without success. 

A possible reason for this failing is suggested by Figure El (left). At fixed- 
Q 2 , the poor \ 2 arises because the fit has much too flat an energy dependence 
at the larger :c-values for all except the lowest Q 2 value. This could be 
corrected by increasing the proportion of the hard term, but this necessarily 
would lead to a steeper dependence at the lower x- values at all Q 2 . This 
interpretation is confirmed by the solid curve in Figure El (left), which shows 
a the result of fitting only to data in the x-range 5 x 10 -4 < x < 10 -2 , and 
then extrapolating the fit to lower x-values, corresponding to higher energies 
at fixed-Q 2 - As can be seen, this leads to a much steeper dependence at these 
lower x- values than is allowed by the data at all Q 2 . An obvious way to solve 
this problem is by introducing saturation at high energies, to dampen this 
rise. 

We note that the aforementioned problem could also be resolved by re- 
stricting the fitted region to x < 10" 3 . In our view, there is no justification 
for this, since the non-diffractive contributions are already small at x ~ 0.01. 
Furthemore, they make a positive contribution which decreases in size as x- 
decreases so that, to the extent that they are not completely negligible, they 
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Figure 2: Comparison of our new dipole model fits to a subset of DIS data 
[T3] . Left: No saturation fits. FS2004 Regge dipole fit (dashed line) and 
(solid line) a fit of the same model to data in the restricted range 5 x 10~ 4 < 
x < 10~ 2 , extrapolated over the whole x-range x < 0.01. Right: Saturation 
fits. FS2004 saturation fit (solid line) and the CGC dipole model (dot-dashed 
line) 



make the problem worse". 



3.2 Saturation fits to DIS. 

Saturation can be introduced into our new model by adopting a device pre- 
viously utilized in [23]. Instead of taking tq to be fixed we now determine it 
to be the value at which the hard component is some fixed fraction of the 
soft component, i.e. 

(r(x m ,r )/a(x m ,r 1 ) = f (13) 

and treat / instead of fitting parameter. This introduces no new 

parameters compared to our previous fit. However, the scale r now moves 
to lower values as x decreases, and the rapid growth of the dipole cross- 
section at a fixed value of r begins to be damped as soon as r becomes 
smaller than r. In this sense we model saturation, albeit crudely, with r the 
saturation radius. 

^For a simple paramaterization of these non-diffractive contributions, sec Donnachic 
and Landshoff |2H| . 



s 



A H 


0.650 


\h 


0.338 


A s 


58.42 


As 


0.0664 


ro 


0.872 


7*1 


4.844 


m 


0.223 







(12) 



Table 1: Parameters for the FS2004 Regge model in the appropriate GeV 
based units. 



The best fit parameter values for what we refer to as our FS2004 satura- 
tion fit are listed in Table El corresponding to a x 2 /data point of 155/156. The 
corresponding fits to the data are shown in Figure |2 (right). Also shown are 
the very similar results obtained using the more sophisticated CGC model 
(no charm fit) j^j. 



A H 


0.836 


X H 


0.324 


A s 


46.24 


A s 


0.0572 


7*1 


4.48 


/ 


0.129 


m 


0.140 







(14) 



Table 2: Parameters for the FS2004 saturation model in the appropriate GeV 
based units. 

It is clear from these results that the introduction of saturation into the 
model immediately removes the tension between the soft and hard compo- 
nents which is so disfavoured by the data. However, it is important to note 
that this conclusion relies on the inclusion of the data in the low Q 2 region: 
both the Regge and saturation models yield satisfactory fits if we restrict to 
Q 2 > 2 GeV 2 , with x 2 /data point values of 78/86 and 68/86 respectively. 

3.3 Predictions for DDIS 

We finally seek support for our suggestion that saturation dynamics may 
already have been observed in the HERA data by turning to a comparison 
to data on diffractive deep inelastic scattering (DDIS). 
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In principle, if evidence for saturation is seen in DIS data, it should 
also be detectable at some level in DDIS data. As we have seen, there is 
a characteristic difference in the predictions of the Regge and saturation 
models for the energy dependence of the DIS structure function F2 at fixed 
Q 2 . If DIS and DDIS are described by the same dipole cross-section, then 
there should be corresponding differences in the predictions for the energy 
dependence of the diffractive structure function F 2 D ^ at fixed Q 2 and fixed 
diffractively produced mass Mx- In the modern parlance, we shall examine 
the x P = (Q 2 + M X )/(Q 2 + s) distribution at fixed (3 = Q 2 /(M X + Q 2 ) and 
fixed Q 2 . 

A recent discussion of DDIS in the context of the dipole model, including 
the predictions of the CGC model, has been presented in reference jH] and we 
refer to that paper for the relevant formulae and more detailed discussion. 

Our predictions involve no adjustment of the dipole cross-sections used 
to describe the F 2 data. However, we are free to adjust the forward slope 
for inclusive diffraction, b, within the range acceptable to experiment and 
in all cases we took b = 4.5 GeV -2 . This simply influences the overall 
normalization of F 2 D ^ which is therefore free to vary slightly. We are also 
somewhat free to vary the value of a s used to define the normalization of the 
qqg component, which is important at low values of (3. In all cases we choose 
a s = 0.1. 

In Figure 01 we show the predictions of our new FS2004 Regge and satu- 
ration models for the xp dependence of the structure function F 2 D ^ at fixed 
Q 2 and (3, together with the corresponding predictions of the CGC model. In 
doing so, we have chosen to focus on (3 values in the intermediate range where 
the predictions are relatively insensitive to the qqg term and to the large r 
behaviour of the photon wavefunction, which are both rather uncertain. 

There is, as expected, a characteristically different energy dependence of 
the Regge model and the two saturation models. There is a hint that the 
data prefer the saturation models, but more accuracy would be needed in 
order to make a more positive statement. 

In Figure E] we show the (3 dependence at fixed xp and Q 2 . Although this 
is unlikely to exhibit saturation effects in a transparent way, it is nonetheless 
a significant test of the dipole models discussed, because different size dipoles 
enter in markedly different relative weightings to the DIS case. For (3 > 0.4, 
where the qq dominates, it is clear that both the Regge dipole and both 
saturation models are compatible with the data, given the uncertainty in 
normalisation associated with the value of the slope parameter alluded to 
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Figure 3: Comparison of the predictions of the FS2004 Regge fit (dashed 
line), the FS2004 saturation fit (solid line) and CGC fit (dot-dashed line) to 
the data on F2 ■ Preliminary data from 
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Figure 4: Comparison of the predictions of the FS2004 Regge fit (dashed 
line), the FS2004 saturation fit (solid line) and CGC fit (dot-dashed line) to 
the data on F^ 3 \ Preliminary data from |14j . 



above. At the lowest (3 values, where the model dependent qqg component 
is dominant, there is a discrepancy between the data and the predictions of 
all three models, which can only be removed by reducing the value of a s well 
below the value assumed " . 



4 Conclusions 

The colour dipole model offers a unified picture of diffractive photoprocesses 
over a wide range of Q 2 . In this paper, we have attempted to throw light on 
the question of saturation in dipole models using both DIS and DDIS data. 

To do this, we have used a particular dipole cross-section in which satura- 
tion effects can be included, or not included, in a simple way. When we speak 
of saturation we refer specifically to the effect whereby the rapid increase in 
energy of the contribution associated with small dipoles is damped at high 
enough energies. If there is no such damping, we say there is no saturation. 

"The relative size of the qq and qqg components is given in more detail in Figure 6 of 



12 



Without saturation, we find that the model is unable to give a satisfactory 
account of the data, although the fit is better than with our original Regge 
dipole model ^U] and many variations on it. Furthermore the discrepancies 
between data and fit are qualitatively of the sort which one might expect 
that saturation effects could remove. On incorporating such effects, we in- 
deed obtain a very good fit to the data, without extra parameters, and the fit 
is of comparable quality to that obtained using the Colour Glass Condensate 
Model [3], which incorporates a more sophisticated treatment of saturation 
effects. Both models yield very similar, successful, predictions for the DDIS 
data in the regions where the qq dipole contribution dominates. 

Finally, we note that the evidence we have presented for gluon saturation 
rests upon the applicability of the dipole model. It is of course well known 
that DGLAP fits to the DIS data work well down to Q 2 ~ 1 GeV 2 [2311)]. 
Such observations are in accord with our observation that the higher Q 2 data 
can be fitted without invoking saturation, and it is only on admitting the 
lower Q 2 data that the evidence emerges. However, Thorne has recently 
pointed out that DGLAP corrections to the dipole formulation are surely 
needed in order to make precise comparison to the data and therefore that 
any claims for saturation are to be viewed with caution [27] . This may well be 
a valid point and it will be interesting to see if our statement that saturation 
is responsible for alleviating the tension between the low and intermediate 
Q 2 data survives a more careful analysis. 
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